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VELOCITY STRUCTURE IN THE ORION NEBULA. I. SPECTRAL MAPPING IN LOW-IONIZATION LINES 
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ABSTRACT 

High-dispersion echelle spectroscopy in optical forbidden lines of 0°, S + , and S 2+ is used to construct 
velocity-resolved images and electron density maps of the inner region of the Orion nebula with a 
resolution of 10 km s^ 1 x 3" x 2". Among the objects and regions revealed in this study are (1) the 
Diffuse Blue Layer: an extended layer of moderately blue-shifted, low-density, low-ionization emission 
in the southeast region of the nebula; (2) the Red Bay: a region to the east of the Trapezium where 
the usual correlation between velocity and ionization potential is very weak, and where the emitting 
layer is very thick; and (3) HH 873: a new redshifted jet to the southwest of the Trapezium. 

Note: A version of this paper with full-resolution figures can be obtained from 
http : //www. if ront . org/wiki/Spectral_Mapping_Paper_I 

Subject headings: H II regions, ISM: Herbig-Haro objects, ISM: individual (Orion Nebula), ISM: jets 
and outflows, techniques: spectroscopy, star formation 
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1. INTRODUCTION 

The Orion Neb ula is the p rime example of a blister- 
type H ii region (Israel 11978( 1 , in which ultraviolet ra- 
diation from the high-mass Trapezium stars heats and 
ionizes the front surface of the molecular cloud OMC- 
1, giving rise to a "champagne flow" of ionized gas 
(|Tenorio-Tagld [19791: lYorkel Il986t IHennev et all 120051) 
away from the cloud. Although the ba sic empirical mode l 
was proposed more than 30 years ago (| Zuckermanl l 1973). 
the nebula continues to present theoretical challenges due 
to the extraordinarily rich phenomenology that is contin- 
ually uncoveredby advances in observational techniques 
(see lO'Delil l200ll and references therein). The Orion 
Nebula presents one of the best opportunities available 
for studying in detail the star formation process in high 
density environments. A rich cluster of more than a thou- 
sand low-mass stars is assoc iated with the ne bula and 
centered on the Trapezium (|Hillenbrandl 119971 ). There 
is evidence that the star form ation rate has accelerated 
over t he past few million years ([Hillenbrand fc Hartmaiml 
1998), with the youngest stars concentrated towards 
the center of the cluster. Indirect evidence from the 
survival of the proply ds in the core of the nebula 
(jHennev fc O'Delll 1999) indicates that the highest mass 
stars are probably even younger (10 4 -10 5 years). Dense 
molecular gas is primarily concentrated in a north-south 
oriented lane (|Johnstone fe B ally 1999), displaced ~ 30" 
to the west of the Trapezium, containing em bedded sites 
at w hich star formation is still ongoing (Smit h et al.l 
120041 ). A thorough understanding of the structure and 
dynamics of the H n region is necessary in order to un- 
derstand the role of the high-mass stars in limiting the 
star formation process by dispersing the molecular cloud 
and in evaporat ing circumstellar mat erial around young 
low- mass stars (jHester fc Descbll2005l ). 
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Previous global studies of the nebula have been 
mainly based on imaging in emis sion line filters 
dHester et al.l 119911: iPogge et all 119921) or spec t ropho - 
tometry (jBaldwin et alJ 119911 : lEsteban et all 1 19981 ) . 
and therefore only give an average over different emit- 
ting regions along the line of sight. High- resolution 
spectroscopic studies show a wealth of detail with 
many ions sh owing complex mu l ti-component kine- 
matic profiles (iCampbell fc Moorel 119181; IWilson et al.l 
19591: iDeharvend 119731: iGoudis et all 119841: lHanell 119871: 

' 19921: lO'Dell fc Wenl 
Massev fc Mea burr 
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IWen fc O'Delll 1199? 
Massev fc Meaburnl 119951: lO'Dell et alJ 120011) . Al- 
though in principle high resolution optical spectroscopy 
can allow the disentanglement of different emission 
regions, in practice progress has been slow due both to 
the complex nature of the object and the incomplete 
nature of the data-sets available. Global spatio-spectral 
studies of the nebula have gener ally been car ri ed out 
with Fabry-Perot in struments (lO'Dell et all [199713 : 
iRosado et all 120011 : Ide La Fuente et alJ 12003 1. the 
velocity resolution of which (> 30 km s _1 ) is not 
sufficient to study any but the highest velocity flows 
in the nebula. 2 R adio recombination line studies (e.g., 
IWilson et al.| [T997) offer the possibility of higher spectral 
resolution but have been confined to hydrogen lines, 
where the thermal Doppler broadening of ~ 20 km s _1 
largely negates this advantage. On the other-hand, 
most multi-li ne high-resolution (~ 5 km s _1 ) echelle 
spectroscopy (jHennev fc O'Delll 1 1999( 1 has only covered 
restricted areas of the nebula . A n otable exception is 
the recent study bv lDoi et all (|2004 hereafter DOH04), 
which covered the entire inner region of the nebula 
in the [N 11] , Ha and [O m] lines with a resolution of 



In this paper, we present a study 
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2 Earlier Fabry-Perot studies (e.g., Dcharveng 1973) were carried 
out at higher resolution, but employed photographic plates, thus 
making quantitative analysis difficult. 
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in the lines of [S n]A6716,6731, [Oi]A6300, [Sm]A6312, 
which is complementary to the DOH04 study. Whereas 
DOH04 concentrated on the high-velocity emission from 
jets and HH objects, this paper emphasizes the more 
subtle variations in the bulk of the nebular gas, which 
typically moves at transsonic velocities of 5-20 km s -1 . 

The four lines that we study cover a range of ion- 
ization stages from low to moderate ionization. These 
correspond to different spatial zones as a result of the 
ionization stratification that is typically found in ion- 
ized nebulae. The [O i] line is emitted principally by the 
partially ionized gas in the ionization front, whereas the 
[S in] line is emitted only by ionized gas in the interior 
of the nebula, where hydrogen is fully ionized, although 
excluding the most highly ionized interior of the nebula, 
where sulfur is triply ionized. The emission of the [S n] 
lines contains contributions from both the fully ionized 
and partially ionized zones and furthermore the differ- 
ence in critical density of the two lines of the doublet al- 
lows the determination of the electron density. Our set of 
lines is of generally lower ionization than that studied in 
DOH04 since the [S n] and [O i] lines extend to more neu- 
tral regions than does [N n] , whereas [S in] is weak from 
the most highly ionized regions, which emit strongly in 
[Oni]. Previous, low spectral resolution studies of the 
electron density in the Orion nebula using these lines 
have suffered from the problem that the derived densi- 
ties are merely an average of the conditions in different 
regions of the nebula, weighted by the [Sn] emissivity. 
However, by kinematically resolving the different com- 
ponents of the line, we are able for the first time to sep- 
arately determine the electron density in different zones 
of the nebula. 

Some preliminar y results from t his st udy have already 
been presented in IHennev et alj (|2005f ) , where the pat- 
terns of global kinematics in the nebula were compared 
with the results of hydrodynamical simulations of cham- 
pagne flows. The current paper contains a much fuller 
description of our dataset and an empirical analysis of 
the complex structure of the spatio- kinematic maps, con- 
centrating on velocities within ±30 km s -1 of the sys- 
temic velocity and presenting evidence for the existence 
of many new nebular features. Our results on the low- 
velocity HH 528 out flow are reported in a separate paper 
( Hennev et al.ll2006f). wherea s another companion paper 
(Henney fe Garcia- Di'az 2006s, referred to hereafter as Pa- 
per II) will present a more detailed study of some of the 
new features we have identified in the nebula, including 
comparison with high-resolution HST imaging. 

In §[2] we describe our observations and the data reduc- 
tion steps that we followed in order to produce flux- and 
velocity-calibrated spectra and velocity channel maps of 
emission and electron density. In §[3] we use these results 
to empirically derive the physical conditions and kine- 
matics of different features in the nebula, many of which 
have not been previously discussed in the literature. 

2. THE OBSERVATIONS AND DATA REDUCTION 

High-resolution spectroscopic observations of the 
Orion Nebula were obtained from two different obser- 
vatories: one set from Kitt Peak National Observa- 
tory (KPNO), and the other set from Observatorio As- 
tronomico Nacional at San Pedro Martir, B.C., Mexico 
(SPM). The first set of spectra was obtained with the 



echelle spectrograph mounted on t he 4 m telescope at 
KPNO and is described in detail in 10 'Dell et all (|2001[ ) 
and DOH04. The 300 x 0.8" slit was oriented north-south 
on the nebula and exposures were taken at a grid of dif- 
ferent declinations (separation « 2"). This dataset, with 
velocity resolution of 8 km s~ , consists of 37 positions 
for [S n] AA6716,6731 A lines in two disjoint regions of the 
nebula: one in the East, centered on RA = 05 h 35 m 2r!5, 
and one in the West, centered on RA = 05 h 35 m ll?70. 

The SPM observations were obtained in 2002 October 
and 2003 January for [S n] (55 pointings) and in Decem- 
ber 2003 for the [Oi] and [Sin] lines (60 pointings) us- 
ing the M ES-SPM instrument (M anchester Echelle Spec- 
trometer, Mcabur n et al.l l2003f ) attached to the 2.1 m 
telescope in its f/7.5 configuration. For the majority of 
the exposures, we used a 150/^m (= 2 arcsec) slit, giving 
a velocity resolution (FWHM) of 12 km s~ , while for 
some exposures we used a 70/^m (= 0.95 arcsec) slit, giv- 
ing a velocity resolution of 6 km s _1 . The 312"-long slit 
was oriented North-South for most of the observations, 
except for one EW-oriented exposure in each line at a 
position 9.6" south of 9 1 Ovi C. We used the star JW499 
(a = 05 h 35 m 15!89, 6 = -05°23 / 49"92: IJones fc Walkerl 
1988) as a reference point for all of our pointings. For 
each position we took two spectral exposures with du- 
ration of 450 s for the 150 /nm slit (with 2x2 pixel 
on-chip binning) , or 900 s for the 70 /Ltm slit (with no on- 
chip binning). Thorium- Argon lamp spectra were taken 
for wavelength calibration between each slit position. In 
order to establish the exact position of the slit in each 
pointing we took direct "slit images" of short duration, 
in which the diffraction grating was replaced by a mir- 
ror, with the slit removed for half of the duration of the 
exposure. This was necessary since the telescope point- 
ing proved insufficiently reliable to guarantee 2" offsets 
between exposures. 3 

Combining the datasets from the two observatories 
gives 92 NS pointings in the [Sn] doublet, spanning an 
interval of 190.6" in RA. For the [Oi] and [Sin] lines, 
which were observed simultaneously due to their prox- 
imity in wavelength, we have only the 60 NS pointings 
from SPM, spanning a similar interval in RA. 

The initial stages of data reduction (bias removal, flat 
fielding, cosmic ray removal) were carried out using stan- 
dard IRAF 4 tasks. This was followed by rectification and 
first-order wavelength calibration of the two-dimensional 
spectra based on the comparison lamp spectra, also using 
standard IRAF tasks. 

After transforming all the spectra to a common he- 
liocentric velocity frame, we carried out a series of fur- 
ther steps in order to produce internally consistent, well- 
calibrated velocity cubes in each of the four emission 
lines. These were implemented using purpose-written 
Fortran and Python routines. 

1. Continuum emission was removed by fitting a 
quadratic function to each row of each two- 

3 This temporary pointing problem at the SPM observatory has 
since been solved. 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under cooperative agreement with 
the National Science foundation 
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Fig. 1. — Sample two-dimensional slit spectra of [Ol]6300A, [S II] 6731 A, and [S III] 6312 A for four representative slit positions. All 
spatial offsets are measured with respect to 1 Ori C: (a) 68.8" W, (6) 22.2" W, (c) 40.4" E, and (d) 88.4" E. The grayscale negative 
images are each scaled to the maximum intensity in each slit and use a "square root" transfer function, which emphasizes fainter features. 
The color, positive-scaled images are composites of the three emission lines for each slit, with [O i] represented by red, [S II] by green, and 
[S III] by blue. 



dimensional spectrum, including only line-free re- 
gions in the fit. 

2. An astrometric solution was found for each of the 
"slit images" mentioned above, using several of the 
Trapezium cluster stars. This allowed us to accu- 
rately determine the slit position of each exposure. 

3. Due to variations in atmospheric transparency be- 
tween exposures, it was necessary to perform a slit- 
to-slit relative intensity calibration. This was done 
by comparison with the spectrum from the hori- 
zontal slit position. 

4. In the same way, the relative slit-to-slit wavelength 
calibration was refined using the horizontal slit 
spectrum. The line centroid velocity as a func- 
tion of position was calculated for each slit and 
small shifts were introduced in the velocity scale 
to force agreement with the horizontal slit at the 
point where they crossed. The shifts required were 
typically less than 2 km s _1 . 



5. For the [Oi] line it was necessary to remove the 
night-sky component, which occurs at km s _1 in 
the geocentric frame. We had chosen the dates of 
the observations such that the heliocentric correc- 
tion means that this component is not blended with 
the main emission from the nebula. Since the pro- 
file of this component was seen to be constant along 
the length of the slit, it was straightforward to re- 
move by fitting an analytic profile to the summed 
spectrum of each slit (a hyper-Gaussian profile was 
used) . The centroid of the fitted profile provided an 
additional check on the accuracy of the wavelength 
calibration for this line. 

Examples of the resultant two-dimensional spectra af- 
ter following all these steps are shown in Figure [TJ Only 
the longer wavelength component of the [S n] doublet is 
shown. In this and all following figures, heliocentric ve- 
locities, Vq, are used 5 and positions are specified as arc- 
second offsets (x, y) with respect to the principal ionizing 



5 The transformation to 
VW = Vp, - 18.1 km s" 1 . 



'Local Standard of Rest" velocities is 
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Fig. 2. — Natural-weighted channel maps of the line cores, color-coded as shown in the top-left panel, showing the coarse velocity changes 
within ±18 km s — 1 of the peak [Sll] velocity of Vq = +22 km s _1 . For each emission line, an identical intensity scale is used for each 
of the three RGB channels, allowing meaningful comparison to be made between the colors in the different maps. The upper-left panel 
indicates the position of each velocity channel on the integrated spectrum of the entire nebula (logarithmic scale) in the three lines. Left 
to right: [Sill] (turquoise), [Sn] (gold), [Ol] (purple). 
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RA offset (arcsec) RA offset (arcsec) 

Fig. 3. — Same as Fig.[2]but for the blue side of the line profiles, corresponding to heliocentric velocities between Vq = —20 km s" 1 and 
+16 km s — . Since the emission line profiles are falling towards the blue in this velocity range, different normalizations are used for the 
images in the 3 velocity bands, as indicated by the horizontal white lines in the top-left panel. 
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Fig. 4. — Same as Fig. [2]but for the red side of the line profiles, corresponding to heliocentric velocities between Vq = +28 km s~ 1 and 
+64 km s _1 . Since the emission line profiles are falling towards the red in this velocity range, different normalizations are used for the 
images in the 3 velocity bands, as indicated by the horizontal white lines in the top-left panel. 
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star of the nebula, 6' 1 Ori C, which has J2000 coordinate s 
a = 05 35 16.46, 6 = -05 23 23.2 (IHog et al.l[200l . 
Transformations between these offsets and the commonly 
used coordinate-based naming scheme of IP 'Dell fc Wenl 
( 1992J) are given in the Appendix. 

The preceding steps have not significantly degraded the 
spatial or velocity resolution of the original data. How- 
ever, in order to combine the individual slit spectra into 
a three-dimensional data cube, it is necessary to degrade 
the resolution to that of the lowest common denomina- 
tor. To that end, we carried out the following steps: 

6. The instrumental linewidths were homogenized by 
smoothing the SPM 70 /zm spectra by 10 km s _1 
and the KPNO spectra by 9 km s _1 to give the 
same instrumental width as for the SPM 150 fj,m 
spectra (12 km s _1 FWHM). 

7. The slit spectra were binned/interpolated 6 onto a 
uniform grid in RA, declination, and heliocentric 
velocity, with voxel size of 0.6" x 0.6" x 4 km s _1 , in 
order to produce the three-dimensional data cube. 
During this stage, we also applied a "de-jittering" 
procedure to each two-dimensional iso- velocity cut 
in order to eliminate unsightly striping due to resid- 
ual inter-slit variations. For each velocity value, a 
mean intensity of extended diffuse emission 7 was 
calculated for each slit, and these were then fil- 
tered across the slits by a 5-point smoothing kernel 
to produce a correction to the intensity of each slit. 

The resulting data cubes for the lines of the three ions 
are shown in Figures to 2] as velocity channel maps. 8 
Each of the images encodes the emission in three adja- 
cent 12 km s _1 velocity channels as blue, green, and red, 
respectively, moving from more negative to more pos- 
itive velocities, as shown on the logarithmically scaled 
line profile of the final panel. We estimate the spatial 
resolution of the resulting images to be 3" x 2" for [S n] 
and 4" x 2" for [Oi] and [Sin]. Despite our best efforts, 
some artefacts remain from the combination of different 
datasets. These are identifiable as bright or dark vertical 
stripes on the images, which are more prominent in the 
line wings, where the intrinsic signal is weaker. They are 
chiefly caused by imperfections in the rectification and 
continuum removal, or in some cases due to reflections 
caused by a bright star falling in the slit (particularly at 
an RA offset of 0). 

Further steps were necessary in order to photometri- 
cally calibrate our data. We used spectrophotometric 
observations from the literature in order to "pin down" 
the absolute line fluxes at different points of the neb- 
ula by comparing with the velocity-integrated intensity 
that we observe at the same point. For [O i l and [Sin], 
we used the observations of IBaldwin et al.l (1991), who 



6 We used an algorithm that reduces to binning in the case where 
the spacing between the original data is finer than the target grid 
spacing and reduces to linear interpolation in the case where it is 
coarser. The algorithm also takes care of interpolating over "bad 
regions" in certain spectra, which were due, for instance, to spec- 
trograph reflections when the slit passes over a very bright star. 

7 Defined as pixels within 3<r of the mean. 

8 The data will be mad e publicly available via the ft p server 
f tp: //orion.phy . vanderbilt . edu/outgoing/ResearchScientist/Ori 



give the absolute flux of the blend of these two lines, to- 
ge ther with the higher spectral resolution observations 
of IBaldwin eTail pOOOh , who give the ratio of the two 
lines. In both cases the observations are for an EW slit, 
positioned to the W of 6' 1 Ori C, which overlaps with our 
own data over a length of ~ 40". Note that, although 
we obtain the [O i] and [S in] lines in the same spectrum, 
we still need to externally calibrate their ratio because 
of the unknown variation in the spectrograph efficiency 
over the wavelength range of the two lines. Since we per- 
formed all the observations with the same spectrograph 
settings, wc have assumed that the relative efficiency is 
the same for a ll slits, and, indeed, comparison with the 
IBaldwin et al.l data confirms this assumption. 

For the [S n] doublet, a greater var iety of spectrophoto- 
metric data exists in the literature fB aldwin et al.lll99ll 
IPogge et al.lll992tlEsteban et al.lll998l ). HoweverTin this 
case, we have the added complication that our data are 
a combination of results from two different instruments 
at two different observatories. Hence the calibrations 
of the SPM-MES and KPNO slits were performed sepa- 
rately. The density-sensitive ratio of the two [Sn] lines 
is of key importance, so we took particular care to verify 
our velocity-integrated results against those of previous 
authors. 

We have also calculated channel maps of electron den- 
sity and of the line ratios [Or] /[Sir] and [Sin] /[Sir], 
which are shown in Figure [5] and [B] In both cases, the 
sum of the two [S n] lines was used in the denominator 
and for the ratios with [O i] and [S in] the [S n] maps 
were first smoothed slightly in order to better match the 
spatial resolutions. 

3. KINEMATIC-IONIZATION STRUCTURE OF THE 
NEBULA 

In this section we provide an empirical description of 
the nebular structure and kinematics as revealed by our 
spectra. We begin with a general overview of the kine- 
matics and follow with more detailed treatments of some 
of the individual objects and spectral features, both those 
that are already known and those that are described for 
the first time in this paper. We have made use of both 
the isovelocity channel maps and the position-velocity 
spectra when identifying features. The channel maps are 
most useful for taking advantage of the image-processing 
capacity of the human brain in order to reveal large-scale, 
spatially coherent features that are not readily apparent 
in the individual slit spectra. However, it is always nec- 
essary to return to the original position- velocity spectra, 
which in many cases are of higher spectral resolution, in 
order to confirm identifications and make precise mea- 
surements. Figure [7] illustrates in cartoon form the loca- 
tion of all the features discussed below. 

When naming new objects, we have tried to follow 
the established conventions of the literature and to em- 
ploy names that evoke only the two-dimensional appear- 
ance of the feature on the plane of the sky. By delib- 
erately avoiding descriptions that presuppose a certain 
model of the three-dimensional structure of the nebula, 
we hope to avoid the problem of the names becoming 
obsolete with advances in our physical understanding. 
Some widely used examples of the existing nomencla- 
ture are rays — narrow, precisely linear features, which 
onmsEQi&e bright or dark, usually aligned with some source 
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Fig. 5. — Blueshifted isovelocity channel maps of [S ll]-derived electron density (logarithmic negative grayscale), [Ol]/[Sn] ratio, and 
[Siii]/[Sii] ratio (linear negative grayscales). Each velocity channel is of width 12 km s _1 . 
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Fig. 6. — Same as Figure[5]but for redshifted velocities. 
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Fig. 7. — Finding charts for features in our channel maps, (a) Blueshifted features assosciated with the Big Arc (dark blue) and HH 
outflows (light blue), together with the newly identified redshifted jet (red). (6) Bright Bars and newly identified compact bars (yellow), 
(c) Dark regions corresponding to low-ionization lanes (olive green), foreground extinction features (mid gray), and the dark arc (light 
gray), (d) Redshifted features seen in highcr-ionization lines (light pink) and lower-ionization lines (magenta), together with reflection 
filaments seen in the far-red (brown), (e) OB stars (filled black circles) and proplyds (open circles). The shading of the proplyds indicates 
their peak [Ol] velocity from 20 km s _1 (blue) to 31 km s — 1 (red). (/) Superposition of all the features. 



of illumination; bars — bright, approximately linear fea- 
tures; lanes — dark, approximately linear features; arcs — 
curved, more or less elongated features, which may be 
dark or bright. 

3.1. Large-scale nebular properties 
3.1.1. The line core 

The general structure of the emission near the systemic 
velocity of the three lines is shown in the 3-color isove- 
locity images of Figure O Note that the centroid ve- 
locity of the molecular gas behind the nebula is around 
Vq = +28 km s _1 , lying between the green and the red 
channels of this figure. The emission is very clumpy and 



filamentary in [O i] , rather less so in [S n] , and much 
smoother and more diffuse in [S in] . Th e principal fea - 
tures of the nebula are well known — see lOoudisl (fl982T ) 
for a compilation of early work on the subject. 

1. The region of peak nebula brightness, lying just to 
the southwest of the Trapezium, centered around 
(25, —25). This is a region of highly inhomogeneous 
brightness, especially in the lower ionization lines, 
with much structure that is only poorly resolved 
in our observations, but that can be appreciated in 
HST images. 

2. The Bright Bar (e.g., lElliott fc Meaburnl ((1971 ). 
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a prominent low-ionization, approximately linear 
emission feature to the southeast, running from 
(-100, -60) to (90, -170). T his is the ionized por - 
tion of the famous Orion Bar ([Werner et al.lfT976h , 
which is a concentration of dense molecular gas, 
offset to the southeast from the ionized emission. 
The region south of the Bright Bar is brighter in 
[S n] than in either of the other lines. 

3. The Dark Bay, a foreground extinction feature, 
which can be appreciated as a general reduction in 
brightness to the northeast, with a southern border 
running from (-100, 0) to (-30, 20). 

As one passes from [O i] , through [S u] , to [S in] , the line 
core becomes more weighted towards the blue overall, as 
can be appreciated both from the integrated line profiles 
(shown in the top-left panel of the figure) and from the 
general hue of each image. 9 However, this trend is not 
universal across the face of the nebula and two broad 
regions of redshifted emission are present in the [S in] 
map (shown in fig. [7R The first of these lies to the east 
of the Trapezium, and we name this region the Red Bay 
since its general shape seems to mimic that of the Dark 
Bay that lies to the northeast of it. It is discussed further 
in § 13.41 below. The second region lies just north of the 
western end of the Bright Bar, and we name it the Red 
Fan due to its appearance in low ionization lines as a set 
of bright plumes that fan away from that portion of the 
Bright Bar. 

There is considerable agreement in the pattern of ve- 
locity changes between the three different lines. In 
each case, the most blueshifted gas is found in a broad 
C-shaped arc to the north, west, and south of the 
Trapezium, which wraps around the Red Bay. Apart 
from the two redshifted regions mentioned above, re- 
duced blueshifts are also seen in the Bright Bar and 
in various other large-scale filamentary structures that 
criss-cross the nebula. In addition to the bright bars 
that have previously been reported in the literature 
(jO'Dell fc Yusef-ZadehlfeOOOl ). we have identified several 
examples of more compact bars towards the core of the 
nebula, which are labelled in Figure [7]?). Furthermore, 
we have also detected a series of dark lanes, which are 
labelled in Figure [TJc. Owing to the narrowness of most 
of these features, a full investigation of their nature re- 
quires a detailed comparison with high-resolution HST 
imaging, which is beyond the scope of the current paper, 
but is carried out in Paper II. 

3.1.2. The blue flank 

Figure [3] shows further isovelocity images, but this time 
for the blueshifted flank of the emission lines. There is an 
overlap with Figure [2] in the sense that the red channel 
of Figure [3] is the same as the blue channel of Figure El 
All of the lines are falling away sharply towards negative 
velocities in this velocity range, so that, if a naturally 
weighted brightness scale (as in Fig. ^ had been used, 
then all the images would be dominated by the red chan- 
nel. Instead, different scales are used for the 3 channels 

9 In all of Figs. [2] to [4] the relative intensity scales of the 3 
velocity channels is identical between the 3 lines, [Ol], [Sir], and 
[S III] , allowing meaningful color comparisons to be made between 
the lines. 



(with maxima as indicated by white horizontal lines in 
the top-left panel) so as to give a more balanced color 
scheme. As a result, care must be taken in interpreting 
the images since the brightest blue features are 10 times 
fainter than the brightest red features. 

The [O i] line is generally very faint in these velocity 
channels, and concentrated in very localized regions. In 
the most negative (blue) channel, one sees only emission 
from high-velocity Herbig-Haro shocks: HH 201 to the 
northwest, HH 202 to the west, HH 269/507 to the south- 
west and HH 203/204 to the southeast. In the least neg- 
ative (red) channel, one sees both low-velocity Herbig- 
Haro objects, such as HH 528 to the south, and the blue 
wings of the bright bars. 

In the [Sn] and [Sin] lines, the most negative channel 
is again dominated by the HH objects, but the other two 
channels show much more extensive emission. In particu- 
lar, the green channel, centered around V@ = —2 km s _1 
is very prominent in the south of the nebula in these 
two lines. Part of this emission is due to the Big Arc 
(jO'Dell et al.lll997aHDoi et al.ll2004h . a broad, elongated 
swathe of intermediate velocity blueshifted emission, run- 
ning from southwest from (—100, 0) to (—50, —60) and 
then more nearly west to (80, —75) (see Fig. 0a). Al- 
though it is visible in [Sn] (faintly) and [Sin], the Big 
Arc is most prominent in higher ionization lines, such as 
[Om], and it is discussed in greater detail in Paper II. 
In addition to the Big Arc, widespread diffuse emission 
is seen in this velocity channel, being especially strong in 
[S n] . The [S n] density sensitive line ratio map (Fig. [5]) 
show this emission to come from low-density gas, which 
we refer to as the Southeast Diffuse Blue Layer and which 
is discussed further in § 13.61 below. 

3.1.3. The red flank 

Figure [4] shows isovelocity images for the redshifted 
flank of the emission lines, with the blue channel here 
being identical to the red channel of Figure O Again, 
unequally weighted brightness scales have been used to 
compensate for the general fall-off in the line profiles 
away from their core. This time, the brightest red fea- 
tures are a factor of 50 times fainter than the brightest 
blue features. 

Even with this help, for [O i] very little emission can 
be discerned above the noise in the redmost channel 
(Vq — +52 to +64 km s" 1 ), except for a faint glow 
in the extreme northwest. In the [Sn] and [Sin] im- 
ages, on the other hand, faint but extensive emission is 
found, concentrated towards the region southwest of the 
Trapezium. In the [S n] image, which has a significantly 
higher signal-to-noise than the other two, the emission 
can be seen to be concentrated into clumps and arcs (see 
Fig. Of), which usually do not coincide with the bright 
filaments seen at bluer velocities. We believe that this 
emission arises from the scattering of nebular emission 
lines by dust particles that lie behind the nebula. See 
further discussion below in § 13.51 

Most of the jets and outflows in the Orion Nebula show 
only blueshifted emission (see iHenney et alJ 120061 ) , but 
we have detected a hitherto unreported redshifted jet in 
our data, which is clearly visible in the [Oi] and [Sn] 
images of Figured] at a position of (—40, —35). This is 
described in detail in § 13.2.21 below. 

At the more moderate redshift of the blue and green 



12 



Garci'a-Dfaz & Henney 



channels (Vq = +28 to +52 km s 1 ), one sees mainly the 
red wings of the flows from the bright bars (see Paper II) . 

3.2. Herbig-Haro outflows 

The classical M42 Her big-Haro objects HH 201-204 
and accompanying jets (|Q'Dell et al.l 11997b! ) are all 
present in our spectra, together with more recently recog- 
nised outflows such as HH 529 (see Fig. |3|). Although 
these flows do appear in the velocity channels that we are 
considering, the majority of their emission is at higher 
blueshifted velocities (see, for example, Fig. which 
shows the HH 203 jet and the HH 203/204 bowshocks). 
We therefore defer detailed discussion of these objects to 
a later paper and merely note their presence here. 

Other Herbig-Haro flows do not show such high ve- 
locities and the entirety of their emission is contained 
within the velocity range considered in this paper. These 
include the large, low velocity HH 528 flow and two ob- 
jects that are described here for the first time: HH 873, 
a newly discovered redshifted jet, and the South Blue 
Arc, which may represent a leading bowshock of the 
HH 203/204 system. 

3.2.1. HH 528 

HH 528 is a large-scale, low-ionizati qn, low-velocity 
outflow in the core of the Orion Nebula (|Q'Dell fe Ballvl 
U999tlBallv et al.ll2000D . Although it is one of the largest 
outflows in the nebula, its existence was recognized only 
relatively recently, owing both to its lack of high- velocity 
emission and the fact that its bowshock, despite be- 
ing very large and bright, is superimposed on the low- 
ionization emission from the Bright Bar. Previous spec- 
troscopic studies have been unable to kinematically re- 
solve the outflow from the systemic nebular emission. 
The multi-line velocity mapping that we present here al- 
lows the emission from the outflow and its bowshock to 
be cleanly separated for the first time from the more 
redshifted nebular gas, showing it to be blueshifted by 
10 ± 3 km s" 1 with respect to the systemic velocity of 
the molecular cloud. Further details of our radial veloc- 
ity me asurements of this o bject are given in a separate 
paper (jHennev et al.ll2006l ). 

3.2.2. HH 873 

This newly detected redsfifted jet is shown in greater 
detail in Figure [3 The top panel shows an HST WFPC2 
image of the region in the [N n] line (f658n filter) at a res- 
olution of ~ 0.1", which is roughly 10 times better than 
the resolution of our spectral observations. The jet is re- 
solved into an irregular, knotty structure, with long axis 
at PA ~ 70°. The jet body is crossed prependicularly 
by several faint arcs along its length, and its eastern end 
is marked by a larger asymmetric arc, reminiscent of a 
bowshock. 

From the morphology of the jet alone, it is not clear 
whether it is travelling in an eastward or westward di- 
rection. If the large arc is indeed a bowshock, then it 
would argue for an eastward motion. On the other hand, 
a prominent knot at the eastern end of the jet (labelled 
East Knot in the figure) bears some resemblance to a 
proplyd cusp, which might represent the exciting source 
of the jet, arguing for a westward motion, although this 
East Knot is not detected as a stellar source in the near 



TABLE 1 

Proper motions of knots in the HH 873 jet 





Position 


PA 


Velocity, Vt 


Knot 


(") 


(°) 


(km s^ 1 ) 


A 


(-48.3, -35.9) 


96 


81 


B 


(-47.3, -36.4) 


88 


68 


C 


(-46.2, -36.7) 


82 


18 


D 


(-44.2, -37.0) 


29 


07 


E 


(-43.6, -37.1) 


92 


20 


F 


(-42.8, -37.2) 


83 


36 


Error 






±10 



TABLE 2 

Radial velocities and electron densities in the HH 873 jet 



Slit 


x Position 
(") 


Velocity, Vq 
(km s _1 ) 


Density, n e 
(cm" 3 ) 


22 


-48.4 


37.2 


3000 


23 


-46.5 


39.4 


2500 


24 


-42.9 


38.0 


3900 


25 


-40.5 


35.4 


6200 


Error 




±1.0 


±500 



infrared ([Hillenbrand fc Carpenter! I2000T ). implying that 
it is not in fact a proplyd. 

In order to resolve this issue, the proper motions of 
individual knots along the jet have been determined 
(O'Dell, priv. comm.) from a pair of f658n i mages 
(8121SE and 5085S1 in Fig. 1 of boi et all [20021 ). sep- 
arated in time by 5.19 years. The results are given in 
Table [T] and are shown graphically in the lower panel of 
Figure [8] These show conclusively that the jet is moving 
towards the east, with the fastest knot having a speed in 
the plane of the sky of ~ 80 km s _1 (the errors in the 
proper motion measurements are of order ±10 km s~ , 
iDoi et aLl 12002). Furthermore, the speed of the knots 
shows a general increase from west to east. 

The contrast of the jet against the background nebula 
is highest in the [O i] line (see Fig. [I]c) and it is also de- 
tected in [S n] , but it is not seen in [S in] . It is also clearly 
visible in the [Nil] 6584 A line observations of DOH04, 
although it is not commented on by those authors. The 
middle panel of Figure [5] shows a color-coded isovelocity 
map of the redshifted side of the [Nil] line from DOH04. 
The jet emission is relatively brighter in the 42 km s~ 
band and therefore shows up as red in the image, against 
the blue/green of the surrounding nebula. 

We have measured radial velocities and electron densi- 
ties in the jet, using 4 of our [S n] slit spectra that cross 
it. The positions and widths of the slits are indicated 
by gray arrows in the bottom panel of Figure [8j and the 
results are given in Table [21 For each slit, the line pro- 
file in a 2" long section across the jet was fitted by 3 
gaussians, which represent the nebula, the jet, and the 
Diffuse Blue Layer. Essentially identical results for the 
radial velocities are obtained from [N n] and [O i] . 

The electron density of the background nebula in the 
same slits ranges from 1500 to 2000 cm~ 3 and it can 
be seen from Table [2] that the jet is significantly denser 
than this, especially at its western end. There is evi- 
dence for an increase in the radial velocity as one passes 
along the jet towards the east, except for the eastern- 
most slit (#22), where the velocity decreases again. The 
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Fig. 8. — Observations of the HH 873 redshifted jet. Top: Portion of the HST WFPC2 image in the [N n] line. Center: Three-channel 
color-coded [Nil] velocity map from ground-based observations, with 12 km s wide channels centered on Vq = 42 (red), 30 (green), and 
18 (blue) km s _1 . Bottom: Finding chart for objects of interest in the image. Also shown are proper motions of individual knots (A-F) in 
the HH 873 jet and the positions of the slits (22—25) used to extract radial velocities. 
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jet component is also weaker by a factor of two in this 
slit, compared with the others, which is consistent with 
the brightness of the jet as seen in the HST image. In 
the westernmost slit (#25), the jet component becomes 
more spatially extended, corresponding with its broader 
appearance on the HST image at that position. No de- 
tectable redshifted component was found in slits that 
cross the East Knot or the bowshock. The jet component 
in all the spectra was found to have an intrinsic width 
(FWHM) of ~ 9 km s _1 after correcting for instrumental 
broadening. 

It is impossible to compare individual slit positions 
with individual proper motion knots due to the wide 
disparity in the spatial resolution. Instead, we com- 
pare only the average radial and tangential velocities 
of the jet. We eliminate proper motion knot D from 
this average because its anomalously low velocity would 
cause its emission to be hidden inside the nebular back- 
ground component in our spectra. The average tangen- 
tial velocity is thus Vt — 51 ± 12 km s _1 in direction 
PA = 89 ± 3° . To calculate the mean radial velocity, it 
is necessary to subtract the velocity of the stellar clus- 
ter, V Q = 28 ± 2 km s^ 1 (jJones fc Walker! 1 19881 ). giving 
Vr = 9.5 ±2 km s _1 . 

Since the jet is not exactly straight, it is hard to say 
with any certainty what its exciting source may be, al- 
though some possible candidates are labelled in the lower 
panel of Figure [8] The dark silhouette proplyd disk 
D183-405 lies close to the western end of the long conti- 
nous section of jet, but the position angle of the ellipti- 
cal disk silhouette is not perpendicular to the jet axis, as 
would be expected on theoretical grounds. Besides, there 
is some indication that the jet extends further to the west 
than D183-405, passing close to the star JW575, which 
is another possible source. However, the jet may extend 
much farther to the west even than this, as evidenced by 
a broken sequence of faint red knots in the middle panel 
of Figure [51 which terminates in the general area of the 
bright proplyd P159-350. This proplyd itself is unlikely 
to be the sou rce, since it is kno wn to have a jet pointing 
to the NNE (jBallv et al J 120001 ). but it also has a fainter 
binary companion which might be the source. Unfortu- 
nately, any attempt to trace the jet over large distances 
is hampered by the presence of many other high speed 
flows in the same general area, such as HH 529. 

3.2.3. South Blue Arc 

Starting roughly 10" south of the nose of the HH 204 
bowshock, we detect a very faint arc of blueshifted [S u] 
emission, which runs about 30" to the south, while bend- 
ing slightly to the west (see bottom- left panel of Fig. [3|) . 
This arc is most visible around —10 km s~ and is dis- 
tinguished from the emission of the diffuse blue layer, 
both by its slightly bluer velocity and by its density of 
~ 1000 cm~ 3 , which is roughly twice that of the diffuse 
blue layer, as can be seen in the top-left panel of Fig. [5]). 

This blueshifted emission arc may possibly represent 
a third bowshock due to the same flow that gives rise 
to HH 203/204. It is probably associated with a co- 
coon of faint, high-ionization emissi on detected in deep 
HST imaging of the same region (|Hennev et alj 120061 
§ 3.5). DOH04 present arguments against the possibility 
of a common origin for HH 203 and HH 204 but these 



TABLE 3 
Proplyds identified in our spectra 





Position 


V peak ([Oi]) 


ID 


(") 


(km s" 1 ) 


P155-338 


(+14, -14) 


+23 


P159-350 


(+07, -26) 


+23 


P167-317 


(-03, +06) 


+30 


P168-326 


(-05, -02) 


+31 


P170-337 


(-08, -13) 


+30 


P177-341 


(-18, -17) 


+25 


P 180-331 


(-22, -07) 


+26 


P182-413 


(-25, -49) 


+25 


P206-446 


(-61, -82) 


+20 



are not conclusive. We suggest that the feature shown 
in Fig. 10 of DOH04, rather than being a second jet 
that drives HH 204, might instead be associated with 
the HH 203 jet. In particular, the knots at 204-448 and 
199-447 seem to form part of one side of the nose of a 
fourth bowshock, with the other side formed by the fea- 
ture that branches up at PA ~ 330° from 204-448. This 
could either be an internal working surface or due to an 
interaction of the jet with the Big Arc. It is worth not- 
ing that the highly blueshifted emission from the HH 203 
jet does show a bright knot at this position (198-441 in 
Fig. 9 of DOH), together with several abrupt changes in 
direction along its length, tending towards smaller posi- 
tion angles as one approaches the source of the jet. The 
progression of radial velocities towards greater blueshifts 
from the South Blue Arc, through HH 204, to HH 203 
suggests that the jet direction is also moving closer to 
the line of sight. Such structures are simila r to those 
seen in simulations of episodic, precessin g jets (|Limll200il : 
iCerqueira fc de Gouveia Dal Pinoll2004f ). 

3.3. Proplyds 

We can identify several of the bright proplyds in our 
data, which are listed in Table GH together with the peak 
[O i] velocity of each proplyd, as determined from our 
spectra. The proplyds are best seen in [O i] since the 
high electron densities (up to 10 6 cm -3 , iHennev et alJ 
2002) seen in these objects imply that the [S n] and [S in] 
lines suffer from collisional deexcitation. They show up 
well in the [Oi]/[Sii] ratio, particularly on the blue side 
of the line. A few of the proplyds are detected in [Sin], 
where they tend to show greater line widths than the 
surrounding nebula, but none are convincingly detected 
in [S n] . The positions and velocities of these proplyds are 
also shown in Figure [TJs. Due to confusion and the effects 
of scattered light, the proplyds close to the Trapezium 
are harder to spot in our data, particularly those to the 
west, and as a result some of the brightest proplyds (e.g., 
L V 5) are missing fr o m ou r sample. 

IHennev fc O'Delll (|1999f ) showed that, although the 
[O i] line can show blue or red wings due to the pho- 
toevaporation flow, the peak [O i] velocity is indica- 
tive of the radial velocity of the proplyd's central star. 
For our sample of proplyds we find a mean velocity of 
Vq — 26 km s _1 and a one-dimensional velocity disper- 
sion (after correcting for an estimated 2 km s~ mea- 
surement error) of a = 2.9 km s~ . This latter is not 
significantly different from the value of 2.5 km s _1 ob- 
tained from proper motion studies of the star cluster 
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pones fc Walker! H988I ). The mean blueshift of the pro- 
plyd stars with respect to the background molecular gas 
(at Vp, = 27-29 km s -1 ) is much less than was found 
by iHennev & O'Delll {l999) for a smaller sample. To 
date, only 6 proplyds have reliably measured inclination 
angles, which req uire fitting kinematic models to multi- 
line spectroscopy (IHennev fc Q'Delll [l999: Hc nnev et all 
2002; iGraham et al.ll2002[ ). By combining our data with 
the [N n] , [O m] , and Ha datasets of DOH04 it should be 
possible to extend this technique to a much larger sample 
and hence determine the three-dimensional distribution 
of proplyds in the nebula. 

3.4. The Red Bay 

This is a large region of relatively redshifted [S in] emis- 
sion, most prominent in the +25 to +35 km s~ range 
(red emission in Fig. [2] and blue emission in Fig. [4]), 
which extends out to the east and east-southeast from 
the Trapezium. On the blue side of the line profile 
(< +20 km s~ ), the Bay region is almost devoid of 
[Sin] emission (see lower-right panel of Fig. [3]), whereas 
on the far-red side it also corresponds to a minimum in 
the brightness of the scattered light (§ 13 . 5|) . All these 
features are illustrated in Figure (7]i. 

In this roughly elliptical, 120 x 80" region, the blueshift 
between the high-ionization lines and low-ionization lines 
is reduced to ~ 2 km s _1 , as compared to the ~ 
10 km s^ 1 that is more typical in the west of the neb- 
ula. The fact that the Red Bay has the appearance of 
an empty cavity in blueshifted velocity channels might 
be taken to imply that the region could be explained 
by the action of the stellar wind from the Trapezium 
stars in confining the ionized gas to a thin layer near 
the ionization front. However, despite its appearance in 
position-velocity space, the Bay is not a thin layer but 
rather has an effective thickness along the line of sight 10 
of ~ 10 18 cm, which is comparable to its lateral extent 
on the sky. 

We therefore propose that the Red Bay is a relatively 
inert region of the nebula, which partakes only weakly 
in the general champagne flow. Maps in molecular lines 
and far-infrared dust emission (jJohnstone fc Bally|[l999l) 
show that there is relatively little dense molecular gas 
behind the nebula at this position, which may have per- 
mitted the ionization front to propagate further along the 
line of sight away from the observer. This region would 
therefore be bounded from behind by a concave, weak-D 
ionization front, the bowl-like shape of which would be 
expected to trap an exten ded zone of slow-mo ving ion- 
ized gas (see, for example. iMellema et al.ll2006[ ). 

3.5. The red scattered component 

In the far-red wing (> +45 km s _1 ) of the lines we 
detect a pattern of bright and dark lanes with a mor- 
phology that is almost identical between the 3 different 
lines. It can be seen as the orange/red emission in Fig- 
ure [4] and is more prominent in [S n] and [S in] than in 
[O i] . This probably represents nebular emission that has 

10 The effective thickness (Hcnnev ct al. 2005) is calculated by 
dividing the emission measure (calculated from the radio free-free 
brightness) by the electron density (measured from our [S II] obser- 
vations). 
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Fig. 9. — North-South profiles of electron density (upper panel) 
and [Sir] surface brightness (lower panel) from a slit situated 40.5" 
E of (^Ori C. The black lines show emission at velocities around 
that of the blueshifted line component (—7 to +7 km s -1 ) while 
the gray lines show emission from the principal line component (+7 
to +50 km s" 1 ). 

been scattered by dust in the PDR and molecular cloud 
behind the ionization front. Because the emitting gas is 
moving away from the dust in the ba ckground cloud, the 
scattered emission will be redshifted dLerov fc Le Borgnd 
fll8^IO'Den"et^[l99l[Hennevlll998l) . FigurelTH shows 
the general layout of the more prominent lanes and 
also indicates (dotted line) the boundary of an extended 
plateau of far-red emission, which is found just to the 
east of the Trapezium. The relationship between the 
scattered component and the bright bars and dark lanes 
is shown in Figure [TJf. 

3.6. The Southeast Diffuse Blue Layer 

A notable property of the [Sn] profiles in the south- 
eastern region of the nebula (Fig. [TJcc?) is that they are 
double-peaked, showing a well-defined secondary emis- 
sion component centered on V© ~ —2 km s , which is 
20-25 km s _1 to the blue of the main component. This 
component is not detected in [O i] or [S in] , although to 
the north of the Bright Bar it merges with emission from 
the Big Arc, which is also visible in [Sin]. The spatial 
distribution of this component, which we will call the 
Diffuse Blue Layer, can be appreciated as the extended 
green emission in the lower left corner of the [S n] image 
of Figure[3] The electron density of this layer is very low, 
as can be seen from Figure [5j 

In the lower panel of Figure [9] we show brightness pro- 
files for red and blue velocity ranges of the [Sn]673lA 
line for a representative slit in the East of the nebula. 
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The blue channel covers velocities —7 to +7 km s~ 
while the red channel covers velocities +7 to +50 km s _1 . 
The Diffuse Blue Layer is the dominant contribution to 
the blue channel brightness south of an offset of —100. 
The upper panel shows that the density of the Layer is 
roughly constant at n c ~ 400 cm~ 3 . The Bright Bar 
crosses the slit for positions between —75 and —100 and 
the brightness of the systemic nebular component be- 
comes so large that its wing also dominates the bright- 
ness in the blue channel. Slightly further to the north, 
a peak in blue channel brightness occurs where the Big 
Arc crosses the slit. Therefore, the density rise in the 
blue channel that is seen between positions of —100 and 
—25 is not due to the Diffuse Blue Layer itself. To the 
north of the Big Arc, between offsets of —25 and 0, the 
Diffuse Blue Layer again dominates the blue channel and 
the density returns to ~ 500 cm -3 . Further to the north, 
the intensity of the Diffuse Blue Layer declines and the 
wing of the systemic component returns to dominate the 
blue channel. 

To the west of the nebula, the Diffuse Blue Layer is 
much less prominent and is only detectable to the south 
of the Bright Bar. A dark band appears to cross the Dif- 
fuse Blue Layer in the —12 km s _1 channel, roughly par- 
allel to the ionization front of the Bright Bar but about 
15" farther away from the Trapezium. 

One of our [S n] slits reaches 40" farther north than 
the area covered by our maps and shows that a similar 
feature to the Diffuse Blue Layer reappears around (+49, 
+ 156). The density is n c ~ 500 cm -3 but the velocity 
is slightly redder at V@ ~ +9 km s~ , although it is still 
well separated from the systemic component, which lies 
at +24 km s _1 with n c ~ 1300 cm" 3 . 

We do not detect this c omponent in [O i] or [S in ] but it 
can also be seen in [Nil] (Henney & O'Dcll 1999) and in 
[On] 3727, 3729 A (|Joneslll992[) . The iHennev fc O'Delll 
observations show it to extend at least another 20" east of 
the region covered by our maps. It seems likely that this 
component is the same as that identifi ed many years ag o 
as a region of line splitting in [Nil] (jDeharvend 119731 ). 
The region ar ound and to the south of the Bright Bar 
coincides with iDeharvenef s Region A, while the similar 
feature found to the North is part of her Region B. A 
third region of splitting (Region C) in the far west of the 
nebula lies outside the area of our observations. 

The lack of [O i] emission suggests that this emission 
does not come from an ionization front, but rather from 
an extended region of fully ionized gas. The lack of [S in] 
emission implies that helium is neutral and that the ion- 
izing spectrum is rather soft. It is therefore possible that 
the Blue Layer is due to a separate region along the line 
of sight that is ionized by # 2 Ori A instead of by 1 Ori C. 
This is also consistent wi th the spatia l distribution of the 
layer in the south east (|Deharvensf s Region A). How- 
ever, the larger s cale distrib u tion o f the [Nil] line split- 
ting mapped by iDeharvengl (|1973| ) argues instead that 
the Blue Layer is instead related to the large-scale cham- 



pagne flow from the nebula. 

4. CONCLUSIONS 

We have mapped the Orion Nebula in a variety of op- 
tical emission lines in order to produce images of the re- 
gion that are resolved in both velocity and in ionization. 
We have also produced the first velocity-resolved elec- 
tron density maps of the nebula. The spatio-kinematic 
information revealed by our maps sheds new light on pre- 
viously known nebular features and also uncovers several 
new components of the nebula. In this paper, we con- 
centrate on low ionization emission within ±40 km s _1 
of the systemic velocity and our main conclusions are: 

1. We find a region to the east and southeast of the 
Trapezium, which we call the Red Bay, where the 
usual correlation between velocity and ionization 
potential is very weak or absent and where the 
emission layer is much thicker along the line of sight 
than in the west of the nebula. This is gas that is 
not sharing in the general champagne flow of the 
nebula, presumably beacause of the local concave 
geometry of the ionization front in the region. 

2. We find an extensive layer of low-density, low- 
ionization, blueshifted emission in the southeast 
of the nebula, which may be ionized by the softer 
spectrum of the nearby 09.5 star # 2 Ori A, or may 
represent the near side of the champagne flow from 
the nebula. 

3. We detect a new redshifted jet, HH 873, to the 
southwest of the Trapezium, and a dense low- 
ionization shell that may represent an outer bow- 
shock in the HH 203/204 flow. 
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APPENDIX 

TRANSLATION BETWEEN THE COORDINATES USED IN THE CURRENT PAPER AND THOSE OF O'DELL & WEN (1994) 

In the literature it is common for objects in the Or ion nebula to be desig nated by a six-digit coordinate-based code 
of the form XXX- YZZ, which was first introduced by 10 'Dell fc Wer] (|1994[ ) and which is based on a truncated form of 
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their J2000 positions, rounded to 0.1 seconds in RA and 1" in declination (e.g., DOH04). For example, in this system, 
1 Ori C would have the designation 165-323. In the current paper, we do not use this system, preferring instead to 
simply give arcsecond offsets (x, y) with respect to 1 Ori C. 

In order to ease the comparison of our results with those of other papers, we here pr ovide equation s to transform 
between the two coordinate systems. The transformation between our (x, y) offsets and lO'Dell fc Wenl 's system is 

XXX = [164.6 - 0.6696 x\ 
Y= [(203.2 -y)/60j 
ZZ= [203.2 -y\ -60Y, 

in which the notation [■ • -J denotes the integer part, padded with leading zeros to the requisite number of digits. The 
inverse transformation is 

x = 245.8 - 1.493 XXX 
y = 203.2 - (60Y+ZZ). 

Note that these equations are only valid for right ascensions between 05 h 35 m and 05 h 36 m and declinations south of 
—05° 20', which corresponds to —650 < x < 246 and y < 203. All our observations fall well within these limits. 
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